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^ : We consider solutions of the four dimensional Einstein- Yang-Mills system with a negative cos- 
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mological constant A = — 3p^, where g is the nonabelian gauge coupling constant. This theory 
corresponds to a consistent truncation of M — 4 gauged supergravity and may be uplifted to d = 11 
supergravity. A systematic study of all known solutions is presented as well as new configurations 
bJQ' corresponding to rotating regular dyons and rotating nonabelian black holes. The thermodynamics 

^ ' of the static black hole solutions is also discussed. The generic EYM solutions present a nonvan- 

, ishing magnetic flux at infinity and should give us information about the structure of a CFT in a 

background SU(2) field. We argue that the existence of these configurations violating the no hair 
conjecture is puzzling from the AdS/CFT point of view. 
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, Recently a tremendous amount of interest has been focused on anti-de Sitter (AdS) spacetime. This 

' interest is mainly motivated by the proposed correspondence between physical effects associated with 

gravitating fields propagatiiig in AdS spacetime and those of a conformal field theory (CFT) on the 
boundary of AdS spacetime [l|, Q . 



' A precise formulation of the AdS/CFT correspondence is made in equating the generating function of 

the correlation functions in the CFT with the string/gravity partition function on the AdS space 0,13 



ZAdsiK *o] = / D [g] D [*] e-^[f.*] ^ ^xp ( f d^x^O^i^.^^]]) = Zcft[K *o]. (1) 

' The induced boundary metric and matter fields are respectively denoted by h and ^'o symbolically, with 

O a conformal operator defined on the boundary of AdSd^i- The integration is over configurations [5, 
of metric and matter fields that approach [ft,, ^0] when one goes from the bulk of AdSd+\ to its boundary. 
. , This conjecture has been verified for several important examples, encouraging the expectation that an 

■ understanding of quantum gravity in an asymptotically AdS spacetime (AAdS) can be obtained by studying 

its holographic CFT dual, defined on the boundary of spacetime at infinity. 

In this context, various AAdS solutions of the Einstein equations coupled with matter fields have been 
studied in the literature. However, although the gauged supergravity AdS theories generically contain 
Yang-Mills (YM) fields, most of the studies in the literature have been restricted to the case of Abelian 
matter content in the bulk. At the same time, a number of results obtained for an asymptotically flat (AF) 
spacetime clearly indicate that a variety of well-known, and rather intuitive, features of self-gravitating 
Maxwell fields are not shared by nonabelian gauge fields. In particular, and in contrast to the Abelian 
situation, self-gravitating YM fields can form particle- like configurations Q. Moreover, the Einstein- Yang- 
Mills (EYM) equations also admit black hole solutions that are not uniquely characterized by their mass, 
angular momentum, and YM charges j^]. Therefore the uniqueness theorem for electrovacuum black hole 
spacetimes ceases to exist for EYM systems. 

As proven by some authors 001 1 even the simple spherically symmetric EYM-SU(2) system with a 
negative cosmological constant A in four spacetime dimensions presents a number of surprising results. A 
variety of well known features of AF self-gravitating nonabelian solutions are not shared by their AAdS 
regular and black hole counterparts. 

Although the picture one finds is very much A-dependent, it is still possible to identify some general 
features. For example, there is always a continuum of regular and black hole solutions in terms of the 
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adjustable shooting parameters that specifies the initial conditions at the origin or at the event horizon, 
rather then discrete points. Depending on the value of A, the spectrum has a finite number of continuous 
branches. Furthermore, there exist nontrivial solutions that are stable against spherically symmetric linear 
perturbations, corresponding to stable configurations with nonavanishing nonabelian charges. The solu- 
tions are classified by nonabelian electric and magnetic charges and the ADM mass. When the parameter 
A approaches zero, an already-existing branch of monopoles and dyon solutions collapses to a single point 
in the moduli space |lfl| |. At the same time new branches of solutions emerge. A general study of these 
configurations together with a stability analysis is presented in Q . 

These spherically symmetric solutions have been generalized in various directions. Axially symmetric 
solutions are discussed in |ni-^S| , NUT-charged configurations and topological black holes with nonabelian 
fields were considered in [ij, [13 . Spherically symmetric, five dimensional solutions of the EYM system 
with negative cosmological constant were examined in llq, this result being extended recently to d > 5 

m 

Although further research is clearly necessary, at least some of the EYM-SU(2) solutions, emerging as 
consistent reduction of d = 11 supergravity on a seventh dimensional sphere [l^ Il9||. have relevance in 
AdS/CFT context. In this case, the ratio between the four-dimensional cosmological constant and the 
gauge coupling constant g is fixed by A/g^ — —3. Apparently the bulk/boundary correspondence for 
AAdS EYM configurations has not received much attention in the literature. Although all such solutions 
containing non-abelian fields are classical they may have a role to play in the full quantum theory. If the 
AdS/CFT correspondence conjecture is indeed correct, it should either be able to account for solutions to 
the EYM system from the CFT viewpoint or be able to demonstrate why the conjecture does not apply 
to them. 

The lack of attention given to AAdS EYM solutions is presumably due to the notorious absence of closed 
form solutions in the presence of nonabelian matter fields in the bulk. (Very few such exact solutions exist 
e.g. [23, which features an effective negative cosmological constant but is not AAdS, and, [2l| for A = 0.) 
However, one can analyse their properties by using a combination of analytical and numerical methods, 
which is enough for most purposes. Euclidean solutions of the A = —3g^ EYM model have been discussed 
in (2^, in a different context, however. The authors of Ref. considered Euclidean wormhole solutions, 
with an conformal infinity. 

The Lorentzian solutions we consider here are very different. Their conformal infinity is the product of 
time and a two-dimensional sphere, plane, or hyperboloid. In the black hole case, these are the nonabelian 
counterparts of the well known AdS4 Einstein-Maxwell solutions. For a, Rx boundary structure, there 
are also globally regular, particle-like solutions. The existence of configurations violating the no-hair 
conjecture is puzzling from the AdS/CFT point of view. Since several distinct solutions with the same set 
of boundary at infinity data may exist, it is not clear how the dual CFT distinguishes between them. 

In the first part of this paper we discuss the features of the AAdS nonabelian solutions with A = 
— 3g^. Although some of these solutions are already known in literature, their properties have not been 
discussed for this particular value of the cosmological constant. New types of EYM configurations with a 
nonvanishing angular momentum are presented as well. All these configurations have a higher dimensional 
interpretation, solving the equations of motion of d = 11 supergravity. 

The second part of this paper attempts a discussion of these solutions in an AdS/CFT context. The 
boundary stress tensor and the associated conserved charges are computed in Section 4, where we also 
present a discussion of black hole thermodynamics. We further argue that studying the A = ~3g'^ EYM- 
AdS system should give us information about the structure of the CFT in a background SU(2) field. We 
conclude with Section 5. The Appendices contain a discussion of technical aspects of solutions' construc- 
tion. 



II. GENERAL FRAMEWORK 

We start by discussing the embedding of the EYM solutions in M-theory. It is known that the standard 
TV = 4 SO (4) gauged supergravity in d = 4 can be viewed as a reduction of d = 11 supergravity on S'^ p^ . 
The bosonic sector of this theory contains two SU(2) fields Ffj_,y, F^^, a dilaton and an axion x with the 
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Lagrangian density 

1 e'^ ~ ~ 1 1 ye'^'^ ~ ~ 

pa panf patiiy papa . , pafiiy papa (n\ 

where the potential V{(l), x) is 

= -25^(4 + 2 cosh (/) + x2e'^). (3) 

It can easily be seen that x = = is a consistent truncation of this theory for F^^i, = F^^u and, as a result 
one finds the action 

\ I d^x^[n -2k- 2Tr{F^,F''-')] (4) 



4 . 

where the field strength tensor = ^^'^Ppiy is 

Fp^ = df,A^ - d^Af, +i[A^, A^] , (5) 

and the gauge field = ^t^A^, with r" an SU(2) basis written in terms of Pauli matrices (the value 
of the gauge coupling g has been set to one without loss of generality). The gauge field transforms as 
A'^ = UAfJJ^ + i{d^jU)U^ under a SU(2) gauge transformations U . 

The value of the cosmological constant as read from V{4) = 0, % = 0) is A = —3. 

As usual, to ensure well-defined Euler-Lagrange field equations, one adds to the action principle Q the 
Hawking-Gibbons surface term 23], /surf = — ^ IdM d^x\/—hK where K is the trace of the extrinsic 
curvature for the boundary dAi and h is the induced metric of the boundary. This term does not affect 
the equations of motion but it is relevant in the discussion of solutions' mass and boundary stress tensor. 

Variation of the action I^J with respect to the metric g'^'^ leads to the Einstein equations 

Rpiy - ^9piyR + ^9pi^ = 21),^, (6) 

where the YM stress-energy tensor is 

T^, = 2Tr(F^„F,;3g"'3 - ^g^.F^pF"''). (7) 

Variation with respect to the gauge field leads to the YM equations 

D^F'^'' = 0, (8) 

where Df_,= df_,+i [A^, •]. 

In this paper we will consider solutions of the EYM equations possesing at least two Killing vectors 
^ — dt, rj — d^, corresponding to a stationary, axially symmetric spacetime. For the time translation 
symmetry, we choose a natural gauge such that the matter fields have no time dependence, dA/dt=0. 
However, a rotation around the z— axis can be compensated by a gauge rotation C^pA = D'i' j2J|, with \1/ 
being a Lie-algebra valued gauge function. This introduces an integer n in the matter ansatz (which is a 
constant of motion) and implies the existence of a potential W with 

Fpp = D^W, (9) 

where W = A^p — ^. The integer n represents the winding number with respect to the azimuthal angle 
if. While if covers the trigonometric circle once, the fields wind n times around. The qualitative features 
of the solutions obtained in the abelian sector of the theory are insensitive to n; however, the winding 
number determines to some extent the properties of the nonabelian solutions. 
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Using the formulas in with the two equal gauge fields, A" = A°', we can uplift any configuration that 
extremizes the action principle Q to d = 11 supergravity. The eleven dimensional metric ansatz reads 

dsli = dsl + + cos^ ^ ^(6'' - A';^dx^' f + sin^ e^Z^®" ~ ^m^^;^)^. (10) 

a a 

The antisymmetric tensor field i^(4) , which appears in the action principle of the d = 1 1 supergravity, can 
be read from [ll[l3 

= -3e(4) + %/2sin^cos^dC A (6 - 9)° A *F^+ 
+ ^ea6c[cos2 C(9 - Af A (9 - Af + sin^ ^(9 - Af A (9 - Af\ A ^F!^ (11) 

where 9", 9° are SU(2) right invariant one forms on two 3-spheres S^^S^. 



III. A = 3 EYM SOLUTIONS 



The field equations ©, ijHIl contain a large variety of solutions. First, any solution (cjp,v, Ap) of the 
Einstein-Maxwell (EM) theory with A = —3 can be viewed as a solution of the EYM equations (|H1) by 
taking A^ = Af_iT where T belongs to the Lie algebra of the nonabelian gauge group |23|. The properties 
of these AAdS abelian solutions have been discussed by various authors (for generic A) , as well as their 
relevance in AdS/CFT context. 

Here we restrict ourselves to the pure nonabelian case. As expected, a number of EYM configurations 
with A = are found to possess AAdS counterparts, with very different properties, however. For A < 0, 
there are also rotating regular EYM solutions, which do not survive in the AF limit. 

In this Section we present a discussion of the EYM solutions which arise as a truncation of A/" = 4 SO (4) 
gauged supergravity, both regular and black hole configurations being considered. Some of these solutions 
have been already presented in the literature, however, without discussing the case of interest A = —3. 
Also, to simplify the general picture we will not consider dyon solutions, except for the rotating case where 
Einstein equations require the presence of a YM electric field. 

The mass-energy of these solutions is computed by using the formalism presented in Section 4. Also, we 
will present here general features of these solutions, without entering into technical details, which are the 
subject of the Appendices A-D. 



A. Regular configurations 



1. The A — —3 Bjoraker-Hosotani monopole solutions 



We start with the simplest case, corresponding to spherically symmetric solutions. These are the AdS 
generalizations of the Bartnik-McKinnon configurations found in 7] by Bjoraker and Hosotani, within a 
metric ansatz 

ds^ = ^+r\de^+sir,Hdip^)-^de (12) 
H{r) p^(r) 

where 

i^(r) = l^^-^r^ 
^ ' r 3 

m[r) corresponding to the local mass-energy density. 

The static, spherically symmetric SU(2) YM ansatz is obtained for a winding number n = 1 and can be 
parametrized by one real function Lo{r) 



A 

2 



i|u)(r)Tid0 + (cot^Ts + w{r)T2)sm9d(f^, (13) 



5 



which impUcs the field strength tensor expression 

F = i|u;'dr A {nde + T2 sm9d(p) - (1 - w^)Tj,de A smOdif^, (14) 

where a prime denotes a derivative with respect to r. For this purely magnetic ansatz, the EYM equations 
take the form 

the equation for p decoupling from the rest, p' jp = — 2w'^/r. 

No exact solutions of this system are known yet, so the equations must be solved numerically. The 
solutions with a regular origin have the following behaviour at r 

w{r) + 0{r^), m{r) ^ 2b^r^ + 0{r^), p{r) po{l - ib^r^) + 0{r'^), (16) 

where b, po are real constants. We are here interested in solutions with AdS asymptotics, which implies 
the following expansion at large r 

Pir) = ^ + ^ + o(^), (17) 



where Wq, M and Ci are constants determined by numerical calculations. M corresponds to the ADM 
mass of the solutions, while determines the value of the magnetic charge, Qm — |1 ~ '^ol (^^^ Appendix 
C). For A = 0, Wo = ±1 are the only allowed values, (M, 6) being restricted to a discrete family indexed 
by the node number of the gauge function ti;(r) 0. 

As found in the AAdS solutions are much less restricted. By varying the parameter b which enters 
the expansion at the origin (|16|) . a continuum of monopole solutions are found. These configurations are 
regular in the entire space. The overall picture depends on the value of the cosmological constant; for 
A = —3, solutions with AAdS asymptotics are found for only one branch with —0.557 < b < 1.31. In 
Figure 1 , we plot the asymptotic quantities M and loq as well as the value of the metric function p at the 
origin as a function of the parameter b. One can see that p(0) diverges as 6 — > hmax (*-e. 5tt(0) — s- 0), while 
M{bmax), ^{bmax) remain finite, while all these parameters appear to diverge as 6 — > b^i^. The critical 
solutions have been studied in Q and are not of interest here. 

A configuration is uniquely characterized by the asymptotic parameters Af, o^o- As a characteristic 
feature of A = —3 globally regular configurations, we notice the existence of zero- and one-node monopole 
solutions only. The nodeless solutions are of particular interest because, as discussed in Q , they are stable 
against linear perturbations. The solution with b ~ 0.619 has wq = and, similar to the 't Hooft-Polyakov 
monopole, has magnetic charge Qm = 1. In this case the AdS boundary conditions play a role similar to 
the Higgs field in the AF case. As seen in Figure 1, for values of b in the interval 0.62 < 6 < 1.31 the gauge 
function uj crosses the r— axis, approaching a negative value at infinity. Although for 0.635 < b < 1.309 
there are two configurations with the same asymptotic value of the gauge function, these are distinguished 
by the value of M, with M{b2) < M{bi) if 62 < fei. 

Typical configurations are displayed in Figure 2 for three diferent values of the parameter b. The solution 
with 6 = 1.3 has = —0.045, corresponding to a near critical configuration with gtt{0) = —0.0032. 



2. Axially symmetric monopoles 



It is well known that, for A = 0, an SU(2) YM-theory coupled to gravity also possesses AF static axially 
symmetric globally regular solutions j26|. labeled by a winding number n > 1. The AAdS generalizations 
of these configurations were reported in featuring very different properties, as expected. Here we 
present a discussion of these solutions for the special case A = —3. 
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As discussed in Appendix A, the minimal axially symmetric YM ansatz is parametrized by four functions 
Hi{T^ 9). These magnetic potentials satisfy a suitable set of boundary conditions at the origin, at infinity 
and on the symmetry axis imposed by finite energy, regularity and symmetry requirements. The spherically 
symmetric YM ansatz is recovered for a unit winding number n = 1, two vanishing gauge potentials 
Hi = H3 = while the other two are equal H2 = H4 = uj{r). 

For large r, the configuration becomes spherically symmetric, with H2 = H4 = wq, the asymptotic value 
of the other two functions being zero. The magnetic charge of these solutions is Qm = ^{1 — ujq\. The 
expression of the gauge ansatz and the boundary conditions are presented in Appendices A, B; see also 
Appendix D for a discussion of the numerical procedure we used to find these solutions. 

The static axially symmetric EYM solutions are obtained within a metric ansatz 

where the metric functions /, m and I are only functions of r and 9. Here r is the radial coordinate, t is a 
global time coordinate, {9, ip) being the usual coordinates on the sphere, with < 9 < n,0 < (p < 2tt. The 
expansion of the metric functions as r — s- cx) is 

fi + f2 sin^ 9 , \ , mi +1712 sin^ 9 , I , h + h sin^ 9 , 1 , , , 

/ = i + + "^ = 1 + ^^:1 + 1^1 + ^^ + o(^),(i9) 

which leads to AAdS solutions. Here /i, /2 are undetermined constants, while li = mi = 2/i/3, I2 = 
6/2/17, TO2 — 14/2/17. As discussed in Section 4, the mass M of these configurations is encoded in the 
parameters /i, /2 

A spherically symmetric spacetime is recovered for / = m and /2 = TO2 = ^2 = (i.e. no angular 
dependence). The coordinate transformation between the resulting line clement and the more familiar 
Schwarzschild-like form H12|l is discussed in Ref. [Tll |. 

The results of the numerical integration indicate that every A = — 3 spherically symmetric regular 
solution appears to present axially symmetric generalizations. For any given winding number, we find 
only one branch of solutions classified by the mass and the value of the parameter loq (we have studied 
solutions up to n — 4). Similar to the n = 1 case, the functions H2 and H4 are nodeless or present one 
node only. The gauge potentials Hi, H3 always present a complicated 0— dependence, while one finds 
usually a small angular dependence for H2, H4. The metric functions /, m and I do not exhibit a strong 
angular dependence, while m and I have a rather similar shape. As expected, the angular dependence of 
the metric and matter functions increases with n. Also, the values at the origin of the metric functions 
decreases with n. The typical profiles of the metric and gauge functions are similar to those exhibited in 
[lT| and we will not present them here. 

In Figure 3 we plot the mass M as a function of the ujq for various winding numbers. We see that the 
n > 1 branches generally follow the picture found for n — I (with higher values of mass, however). 

The energy density of the matter fields e, shows a pronounced peak along the p-axis and decreases 
monotonically along the z-axis (with z — rcos9, p = rsa\9). The contours of equal energy density 
e = — T(* are two-torii and squashed two-spheres. The peak of the energy density along the p-axis slightly 
shifts outward with increasing n and increases in height. For a fixed value of wq, the mass of the solutions, 
M(n), increases with n. For example, for ojq = 0.05, one has M(l) = 1.014, M(2) = 2.467, M(3) = 4.747 
and so on. As a general feature, the particle-like nonabelian solitons are less massive than the extremal 
RNAdS solutions with the same magnetic charge. 

We have found it difficult to obtain axially symmetric generalizations of the spherically symmetric 
solutions near limits of the 6-interval, with large errors for the functions. A different metric paramctrization 
appears to be necessary. 
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3. Rotating regular solutions 



An interesting physical question is whether these static nonabehan regular solutions can be generalized 
to include an angular momentum. For the AF case, contrary to results from perturbation theory p^ . 
no rotating generalizations of the BK solutions seem to exist. In this case, the At components of the 
gauge field act like an isotriplet Higgs field with negative metric, and by themselves would cause the other 
components of the gauge field to oscillate as r — > oo psj, which would imply an infinite mass. Rotating 
solutions are found by including in the theory a triplet Higgs field with a vacuum expectation value greater 
than the asymptotic value of At poj . 

For A < 0, there are no boundary conditions to exclude a nonabelian solution with nonzero electric 
potential; dyon EYM solutions already exist in the spherically symmetric case 0- The existence of dyon 
solutions without a Higgs field is a feature for AdS spacetime; if A > the electric part of the gauge fields 
is forbidden |3,|33|- This makes possible the existence of rota ting regular configurations, too. A discussion 
of these solutions in a fixed AdS background is presented in [13 , where is argued that they survive in the 
presence of gravity. However, the general picture appears to be very complicated, crucially depending on 
the value of the cosmological constant. Here we analyse the solutions' properties for A = —3. 

These rotating EYM solutions are found for a metric form generalizing (|18|l for an extradiagonal metric 
component g^f, which satisfies also the circularity condition [3l| 

ds' = t(7^ + ^'^^') + ^(^<^ + -^^)' - - ^^')'^^'' (21) 

/ -L ^ -^r J r o 

where the four metric functions /, m, / and fl depend only on the coordinates r and 6. At infinity, the 
asymptotic form of the metric functions /, I, m is still given by H19|) while the metric function associated 
with rotation decays as 

+ (22) 

ill 32 being two real constants. As discussed in Section 4, the mass-energy M of the rotating solutions is 
still given by (|2()(l . while their angular momentum is 

The YM ansatz contains in this case six functions: four magnetic potentials Hi [i — 1,..,4) and two 
electric potentials i/s, H^. The boundary conditions satisfied by the magnetic potentials are similar to 
static case. For the electric potentials we impose at infinity = V cosO^ Hq = Vs'mO. The parameter V 
corresponds to the asymptotic magnitude of the electric potential At (oo) and determines the properties of 
these solutions. In the abelian case, by using a suitable gauge transformation one can set set V = (or any 
other value) without any loss of generality. In the nonabelian theory, however, such a gauge transformation 
would render the whole configuration time-dependent. In this case, V enters the asymptotic expansion for 
the nonabelian field strength and has a physical relevance (although it does not contribute to the electric 
YM charge). As discussed in Appendix C, the value V ^ implies a purely magnetic, static configuration 
At = Q = 0. Technical details on these solutions, including the boundary conditions are presented in 
Appendix B, D. 

Globally regular rotating solutions are found by staring with a purely magnetic EYM configuration with 
a given loq and increasing the value of V. Here we consider only rotating solutions with the lowest winding 
number n — 1, although we obtained a number of configurations with n = 2 also. 

The branch structure of the spherically symmetric solutions is preserved in the presence of rotation. 
The rotating configurations depend on two continuous parameters: the value luq of the magnetic potentials 
H2, H4 at infinity (which fixed the magnetic charge Qm = ~ ^ol) magnitude of the electric 

potential at infinity V. 

For a given value of the magnetic charge, a branch of rotating dyon solutions emerges smoothly from 
every spherically monopole solution and extends up to some maximal value of V beyond which gravity 
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becomes too strong for regular dyons to persist. We notice that the value at the origin of the metric 
function / decreases with increasing V and tends to zero as V approaches the critical value Vmax (which is 
Wo dependent), corresponding to the formation of a horizon. With increasing V, the dyon becomes more 
and more deformed. The mass, angular momentum and electric charge increase with V and we find again 
a maximal value for the magnitude of the electric potential at infinity Vmax- Alternatively, we may keep 
fixed the magnitude of the electric potential at infinity and vary the parameter luq. Again, it has proven 
difficult to obtain rotating generalizations of the spherically symmetric solutions near the limits of the 
6-interval. In Figure 4 we present the properties of typical branches of solutions for a fixed value of ujq 
(Figure 4a) and for a fixed V (Figure 4b). 

All solutions we have found present nonvanishing nonabelian electric and magnetic charges, representing 
rotating dyons. A vanishing Qg implies a nonrotating, purely magnetic configuration. However, we find 
dyon solutions with vanishing total angular momentum (e.g. J = for ujq = 0.895, V = 2) that are not 
static (locally T^^O). 

Similar to the static case, the functions H2 and H4 are nodeless or present one node only, although they 
have a small 9 dependence, while Hi, H3 and the electric potentials depend on 6* in a complicated way. 
The metric functions f,l,rn present a rather small angular dependence. 

For all configurations, the energy density of the solutions has a strong peak along the p axis, and 
it decreases monotonically along the symmetry axis, without being possible to clearly distinguish any 
individual component. 



The spherically symmetric black holes found by Winstanley in were the first nonabelian solutions 
with AdS asymptotics presented in the literature. More details on these configurations have been presented 
in 0, including dyonic black holes. These solutions obviously violate the no-hair conjecture and present 
static axially symmetric generalizations that are absent in the abelian sector. 

Also, it is well known that for A < 0, the EM theory has black hole solutions for which the topology of 
the horizon is an arbitrary genus Riemann surface. The EYM-SU(2) counterparts of these solutions with 
a nonspherical event horizon topology have been discussed in . 

The properties of the AAdS nonabelian solutions are strikingly different from those valid in the AF case. 
For example, black holes with A < exist for continuous intervals of the parameter space (the value of the 
gauge field on the event horizon), rather than discrete points. Also, there are configurations for which the 
gauge field has no zeros. Moreover, some of these configurations are stable within a perturbation theory 
approach. 

Since the features of these configurations depend on the value of the cosmological constant, we present in 
this Section an analysis of their properties for A = —3. Apart from some classes of configurations already 
known in the literature, we present here numerical arguments for the existence of rotating black holes with 
nonabelian hair, generalizing for an SU(2) field the Kerr-Newman-AdS solution. 



We start by discussing the better known spherically symmetric solutions and their topological black hole 
counterparts. 

These solutions are obtained within a metric ansatz generalizing H12|l for a nonspherically symmetric 
topology of the event horizon 



B. Black hole configurations 



1. n = 1 static solutions 




H{r) 



dt^ 



(24) 



H{r) 




where 



H{r) = k 



2m{r) 



Ar^ 



(25) 



r 



3 
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Here dfll ~ dO^ + f'^{9)dif'^ is the metric on a two-dimensional surface S of constant curvature 2k. r is 
the radial coordinate for which r ^ oo defines the asymptotic region. The discrete parameter k takes the 
values 1, and —1 and implies the form of the fmiction /fe(0) 

{sind, for k = 1 
e, for k = (26) 
sinh0, for fc = — 1. 

The topology of a constant {t,r) slice is Hg. When k = +1, the universe takes on the familiar spherically 
symmetric form, and the (0, if) sector has constant positive curvature. When fc = 0, the S is a flat surface, 
while for fc = — 1, the {9, ip) sector is a space with constant negative curvature, also known as a hyperbolic 
plane. When S is closed, we denote its area by Vk- Also, we set Vfc/47r = 1 in all numerical data we present 
in this paper. 

The construction of the SU(2) connection is presented in Taking into account the symmetries of 
the line element (|24|) we find 

A = ^{u;{r)nd9 + (^^rg + uj{r)T,) f^d^} , (27) 

which reduces to H13|l for fc = 1. 

As a result, we obtain a simplified YM curvature 



ij^'rifir Ad0 + fkUi'T^dr A dif + (w^ - k) fkT3.de A dip^. (28) 
The EYM equations reduce to 



^ 2 



' '^^ . {^^-kY Hlo' lo{lu^ - fc) p' 2 ,2 

m =uj i?H — , ( ) = , — = — uj , (29) 

Zr^ p pr'' p r 

We find the following expansion near the event horizon which is located at r — > 0, 

m{r) = Y (^k- +m'{rh){r~rh), uir) ^ iUh + lj' {rk)ir ~ r,,), (30) 

where 

, (w^ - fc)2 , rhUJh{ujl-k) 

(since the equations (|29|) are invariant under the transformation lu —>■ —lu, it is enough to consider only 
values of lu^ > 0). For fc = 1, w{r) — ±1 corresponds to vacuum Schwarzschild-AdS solution, while 
u'(r) = is the abelian Reissner-Nordstrom-AdS (RNAdS) solution, with unit magnetic charge. Note 
that, apart from embedded abelian configurations, no extremal solutions with reasonable asymptotics 
exist in this case. 

The condition for a regular event horizon is H' {rh) > and places a bound on LOh 

2m' [r^) = M^:^ < k - Arl (32) 
and implies positiveness of the quantity uj'[rii). For A = —3, this relation implies ujh < ^h(max)i with 

^h{max) =k + rh\J'irl+k, (33) 

(with a minimal event horizon radius rh — I/a/S for fc = —1). The asymptotics as r ^ 00 are still given 
by (|17|l for any value of fc. As argued in Section 4, the configurations' mass is given by MVfc/47r, with M 
the asymptotic value of the metric function m{r). 
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The Hawking temperature of the solutions is evaluated from the surface gravity k as given by Th = 
k/ (27r), with 

(34) 

r=rh 

This implies 

4:TTrhp{rh) 

For every considered value of r^, we find regular black hole solutions with Wh taking values in only one 
interval < lj/i < cj^, where Luf^ is always smaller than LiJh{max)- There are also solutions for which loq > \ 
although LOh < 1. The behavior of the metric functions m and p are qualitatively similar for any value 
of k. However the gauge field behaviour depends on the topology of the event horizon. For k ~ 0,-1, 
in contrast to the spherically symmetric case, we find only nodeless solutions. This can be analytically 
proven by integrating the equation for w, (Haj'/p)' = (^(w^ — k)/pr'^ between rh and r; thus we obtain 
Lu' > for every r > r^- For fc = 1, solutions where lu crosses the axis can exist for small enough values of 
rh- For large values of r/j (e.g. = 1), only nodeless solutions are found in this case too. In this case, for 
sufficiently small ujh, all field variables remain close to their values for the abelian configuration with the 
same r^. Significant differences occur for large enough values of tOh and the effect of the nonabelian field 
on the geometry becomes more and more pronounced. 

In contrast to the picture found in jl5j |. the k = —1, A = —3 configurations always have to > 0. 
The black holes therefore only occur with positive values of the mass. Typical solutions in this case are 
presented in Figure 5. In Figure 6, M, ujq, the value of the metric function p(r) on the event horizon and 
the Hawking temperature are plotted as a function of Uh for fc = 0, ±1 black holes and several values 
of r/j. Note that as the temperature approaches zero {i.e. as the solutions approach extremality) , these 
physical quantities all diverge. 

In Figures 7-8 we plot the mass M and the value of the gauge potential at the event horizon Wh for 
several value of the asymptotic value of the magnetic gauge potential as a function of the event horizon 
radius. For fc = 1, the corresponding solution with a regular origin is approached as 0. For topological 
black holes, we noticed the existence of a minimal event horizon radius Tc, for any given wq- The Hawking 
temperature vanishes as rh—* and a naked singularity develops, while the mass stays finite. 

The discussion in 0| , ^3 on the stability of these black hole solutions within the perturbation theory 
can easily be applied to A = —3. It follows that all fc = solutions are stable; the k — —1 solutions with 
Wo > 1 are also stable as well as the nodeless spherically symmetric solutions. 



«2 = -(l/4).g*V^(5,gtt)(5,ff«) 



2. Static, axially symmetric black holes 

Similar to the regular case, the fc = 1, n = 1 solutions discussed above admit static axially symmetric 
generalizations. (Static fc = 0, — 1 topological black holes with a winding number n > 1 are also likely to 
exist but the corresponding EYM ansatz has not yet been considered in the literature.) 

The situation for a nonabelian field is very different from the EM theory, where the static black hole 
solution is spherically symmetric (or, for A < belongs to one of the three cases (|26|l . with the same 
amount of symmetry) . 

The properties of the AAdS axially symmetric EYM black holes were addressed in j however without 
consideration of the case A = — 3. The metric ansatz in this case is given again by H18|l . with a gauge 
potential Ai written in terms of four functions Hi. We require the horizon of the black hole to reside at 
a surface of constant radial coordinate r — r^, where gtt{fh) = 0. Similar to the procedure in the regular 
case, axially symmetric solutions are obtained by extending the n — 1 configurations to higher values of 
the winding number. These solutions are AAdS and have a regular event horizon but for n > 1 they are 
not spherically symmetric and the event horizon gets deformed away from spherical symmetry. 

The boundary conditions at infinity and on the symmetry axis are similar to those used in the regular 
case; in particular, the asymptotic expansion (|19|l is still valid (see also Appendix B). 
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The surface gravity k turns out to be constant at the horizon, as required by the zeroth law of black 
hole thermodynamics. To derive its expression we use the asymptotic expansion near the event horizon in 

5 = {r - rh)/rh 

fir, e) = f2ie)S^ + 0(5)3, ^(^^ ^ m2{e)5^ + 0(<5)^ l{r, 9) = l2{e)5'' + 0{5f. 
From Th — K/{2n), we finds for the Hawking temperature 

. f^i'^i'^^y^^ , (36) 

which is constant as a consequence of the (r, 0) Einstein equation implying /2TO2,9 — 2m2f2,s- 
For the line element H24|l . the area A of the event horizon is given by 



A = 2.l dOsinO^h^^rl (37) 







which allows a computation of the black hole entropy S — A/A. 

The axially symmetric black hole solutions depend on two continuous parameters (r^, ujq) as well as the 
winding number n > 1. 

The behavior of the solutions is in many ways similar to that of the axially symmetric solitons. Again, 
starting from a spherically symmetric black hole configuration with given wq we obtain higher winding 
number generalizations. Axially symmetric generalizations seem to exist for every spherically symmetric 
black hole solution. For a fixed winding number, the solutions form a branch, which can be indexed by 
the mass and the nonabelian magnetic charge Qm, = "-11 ~ ^ol- This branch follows the picture found for 
n = 1 (with higher values of mass, however). This is in sharp contrast to the A = case, where only a 
discrete set of solutions is found 32]. Also, the Kretschman scalar K — RijkiR^-'^^ remains finite for every 
{r > rh, 0). One finds that the deviation from spherical symmetry increases with growing n. 

Once we have a solution, the horizon variables such as ThtA are calculated in a straightforward way 
from H3t)|) . H37|) . The mass of the solution is computed by using the relation H2U|) . extracting the values of 
the coefficients /i, /2 from the asymptotics of the metric functions. In Figure 9 we plot the mass M, the 
Hawking temperature and the entropy as a function of for black hole monopole solutions with — I 
and n = 1, 2, 3. 

The gauge functions H2, H3, H4 start always at (angle dependent) nonzero values on the event horizon. 
For r/i = 1, we find only solutions where H2, H4 do not cross the r axis. The gauge function H2 is always 
almost spherically symmetric, while the gauge functions Hi and H3 are much smaller than the functions 
H2 and 

For the considered solutions, the metric functions m, /, I do not exhibit a strong angular dependence. 
These functions start with a zero value on the event horizon and approach rapidly the asymptotic values. 
The functions m and / have a rather similar shape, while the ratio m/l indicating the deviation from 
spherical symmetry is typically close to one, except in a region near the horizon. The typical profiles of 
the metric and gauge functions we find for A = —3 are similar to those presented in for other values 
of the cosmological constant. 

The horizon has topology, but geometrically is not a sphere, since its circumference along the equator 
Le turns out to be different from that along a meridian Lp 

Le^ [ d(pj^rsm9 = 27rr,i J , (38) 




r—rh,^— const. 




(39) 



However, for these static solutions, one finds a small deviation from spherical symmetry (as measured by 
the ratio Lg/Lp), at the level of few percent. 
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3. Rotating black holes 

AF rotating hairy black holes in EYM theory were obtained in [s^], within the standard Lewis- 
Papapetrou parametrization of the metric and a gauge field ansatz consistent with the circularity and 
Froebenius conditions. These solutions possess three global charges: mass, angular momentum and non- 
abelian electric charge. Although they possess nontrivial magnetic gauge fields outside the event horizon, 
they do not carry a nonabelian magnetic charge. 

Obviously, the static AAdS solutions should also possess rotating counterparts, representing nonabelian 
generalizations of the Kerr-Newman-AdS solution. However, the construction of such hairy rotating solu- 
tions represents a very difficult task since it involves the solution of a large number of coupled nonlinear 
partial differential equations for the metric and gauge field functions and a much richer set of possible 
boundary conditions as compared to the AF case. 

Here we present the first set of rotating black hole configurations with A = — 3. The ansatz we used in 
this case is similar to that employed for the regular rotating solitons, with the same asymptotic expansion 
as r — > oo. In particular, the expressions H19|) . H22|) are valid in this case too, as well as the expressions (|2Uf) . 
(|23|) of the mass energy and angular momentum. The boundary conditions and details on the numerical 
integration are presented in Appendices B, D. 

For a given winding number, the rotating nonabelian black hole solutions depend on four continuous 
parameters: two geometric parameters - (r/j, O^) representing the event horizon radius and the value of 
the metric function 17 at the horizon respectively, and two parameters associated with the gauge field: 
which gives the magnetic charge and V which is the asymptotic value of the electric nonabelian potential 
At(oo). Not surprisingly, different from the AF case, rotating black holes are found also for V ^Q. 

The complete classification of the solutions in the space of these four physical parameters is a considerable 
task, whose scope is beyond the aim of this paper. We have studied mainly rotating configurations with 
Th — 1 and several values of wq, although a number of solutions have been found for other values of r/j. 
Although rotating black hole solutions should exist for any value of n, we restrict here to a unit value of 
the winding number. 

The properties of the horizon can be computed similar to the static case. The surface gravity is obtained 
from 

^~l/A{D^Xu){D'^xn (40) 

where the Killing vector x — ^~ {^h/rh)v — dt, r] = d^) is orthogonal to and null on the horizon. It 
can be proven that the expansion near the event horizon remains valid in the rotating case, which 
implies the expression for the Hawking temperature. We further consider the area A of the black hole 
horizon computed according to (|37|l . and the deformation of the horizon, quantified by the ratio L^./ Lp of 
the circumferences along the equator and the poles. 

To construct a rotating solution, we start from the corresponding static black hole configuration with 
=0, V = Q. As we increase (fi/i, V) from zero via the boundary conditions, while keeping rh fixed, a 
first branch of solutions forms. For a given V , this branch ends at a critical value Vlh, which depends on the 
value of (o^o, J'/i), and the numerical errors increase dramatically for fJ^ > r2^(cr) rendering the solutions 
increasingly less reliable. As — *■ ^h(cr)^ the geometry remains regular with no event horizon appearing 
for r > r/i, and, the mass and angular momentum approach finite values. As found in 33] for A = 0, 
a second branch of solutions bends backward toward = 0; there the mass and angular momentum 
diverge with $7^^ in the limit Vl^ — > 0. Therefore we expect a similar picture in the AAdS case. However 
the numerical construction of such configurations presents a considerable numerical challenge beyond the 
scope of the present work. Also, the existence of other branches of AAdS rotating solutions, not necessarily 
connected to the static configurations, might be possible. 

The picture gets simpler if we study the dependence of the solutions as a function of V (the magnitude 
of the electric YM potential at infinity) for fixed (cjo, rh,^h)- In this case the solutions share a number 
of common properties with the rotating regular counterparts. In Figure 10 we plot the mass, angular 
momentum, electric charge and the contribution of the electric field to the total mass as a function of V 
for a fixed value of loq. As seen in this picture, the mass, angular momentum and electric charges increase 
with V and we find again a maximal value for the magnitude of the electric potential at infinity. 
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The functions H2 and are always nodeless, although they have a small 6 dependence, while iJi, 
and the electric potentials depend on 0-angle in a complicated way. The metric functions f,l,m present 
a rather small angular dependence, the metric function U, presenting a strong dependence of 9 for small 
values of V. 

For all configurations, the energy density of the solutions has a strong peak along the p axis, and it 
decreases monotonically along the symmetry axis. In contrast to the rotating regular case, we found no 
locally rotating solutions with vanishing total angular momentum, although such configurations are likely 
to exist. 

Further details on these rotating black hole solutions, as well as a discussion of the dependence of the 
solutions on the value of the cosmological constant, will be presented elsewhere. 



IV. A COMPUTATION OF PHYSICAL QUANTITIES 



A. The counterterm formalism 



The mass, angular momentum and action of the solutions discussed in Section 3 is found by using the 
counterterm formalism proposed by Balasubramanian and Kraus 34] to compute conserved quantities for a 
spacetime with a negative cosmological constant. This technique was inspired by AdS/CFT correspondence 
and consists in adding suitable counterterms let to the action. These counterterms are built up with 
curvature invariants of a boundary dM. (which is sent to infinity after the integration) and thus obviously 
they do not alter the bulk equations of motion. 

The following counterterms are sufficient to cancel divergences in four dimensions [35l | , for vacuum 
solutions with a negative cosmological constant (to agree with the standard conventions in literature, we 
set the usual factors I/IGttG in the action principle Q and I/SttG for the Gibbons-Hawking boundary 
term) 



(41) 



SttG 





2 










IdM 







Here R is the Ricci scalar for the boundary metric ft., while P = — 3/A = 1. 

Using these counterterms one can construct a divergence-free stress tensor from the total action 
I=Ihuik+Isurf+Ict by defining 



T — 



SI 



1 



SttG 



Kh,,u - 



(42) 



where Eat is the Einstein tensor of the intrinsic metric hat- The efficiency of this approach has been 
demonstrated in a broad range of examples, the counterterm subtraction method being developed for its 
own interest and applications. If there are matter fields on M additional counterterms may be needed to 
regulate the action (see e.g. js^ for such an example in EYM-dilaton theory). However, we find that for 
a pure SU(2) nonabelian matter content in four dimensions, the prescription H41|) removes all divergences 
(a different situation is found for the five dimensional AAdS nonabelian solutions where the counterterm 
method fails and logarithmic divergences are presented in the total action and the expression of mass 0|). 

Having obtained the boundary energy-momentum tensor, one can determine the conserved charges 
corresponding to the Killing vectors as explained in (34> . i35i |. The usual prescription is to first pick a 
spacelike surface S on the boundary with metric (Tat- The boundary metric is written in the following 
form 



hf.^dx'^dx'' = -N^dt^ + (Tabidx" + N^dt){dx'' + N°dt). 
The conserved charge associated to a symmetry generated by the Killing vector is 



(43) 



(44) 
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where n'^ is a timelike unit normal to S. 

The conserved charge associated with time translation is the mass, the angular momentum being the 
charge associated with the Killing vector d^. 

B. n = 1 static configurations 

We consider first the case of n = 1 spherically symmetric and topological black hole static configurations. 
The results we find by using the asymptotic expressions H17|l for the boundary stress tensor at large r are 

T I 1 £{e^ - Ajwl - m (I 

''^ 8ttG r 327rG ~^ {r^ 



" 8ttG ir 32nG 

Thus the leading order terms in this expression boundary stress tensor are similar to the (topological-) 
Schwarzschild-AdS4 black holes. The presence of nonabelian matter is manifest in the second order of this 
expansion only. 

It can easily be verified that the mass of these solutions computed from (|44|l is given by AjfVfe/47rG, 
where Vk is the area of the surface E. Obviously, these solutions have a vanishing angular momentum. 

C. Axially symmetric configurations 

We consider now the general case of a axially symmetric, rotating spacetime described by the line element 
(I21f) . All relevant expressions can easily be derived by using the the asymptotic form of the metric functions 
(tl^ . if^ . The boundary metric in this case is the Einstein imiverse, ds'^ = £'^{d9'^ + sin^ Odip"^) — dt^ . 

One finds in this way the large-r expansion 

1 /2/i 6/2 . , ^ / 1 



SttG^ r V 3 17 



sm^ ' 



sin^e /2/i , 18/2 . 2.\ , 



T. = -^r^ + ^sin^.Uo^l 



87rG^3 r V 3 17 
^ ((,,+, 2 sin^^)sin^.)+G(-L 



Direct computation shows that this stress tensor is traceless. This result is expected from the AdS/CFT 
correspondence, since even dimensional bulk theories with A < are dual to odd dimensional CFTs that 
have a vanishing trace anomaly. The corresponding expressions for the static case are found by taking 
ji = j2 = in the above relations. For spherically symmetric configuration found within the metric ansatz 
l|TH|) . f2 — and the angular dependence vanishes, as expected. 

By using this relation, we find that the mass and angular momenta of the solutions are given by the 
relations (|20|l and H23|l respectively. For static solutions, the parameter /2 describes the deviation of the 
solutions from the spherically symmetry. 

D. Euclidean action and entropy 



One can use the counterterm expression (|41|l to compute the regularized gravitational action and the 
prove that the entropy of the AAdS hairy black holes is one quarter of the event horizon area. 
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Here we start by constructing the path integral 

Z = J D[g]D[^]e-'^^3M^ (47) 

integrating over aU metrics and matter fields between some given initial and final hypersurfaces, 4* corre- 
sponding to the SU(2) potentials. By analytically continuing the time coordinate t it, the path integral 
formally converges, and in the leading order one obtains 

Z ~ e-'"' (48) 

where Id is the classical action evaluated on the equations of motion of the gravity/matter system. Since 
the Euclidean approach becomes problematic for nonabelian solutions with an electric potential js^ . we 
restrict here to compute the action of static, purely magnetic solutions. We should also remark that the 
variation of the action Q gives the correct equations of motion only if the gauge potential is held fixed 
on the boundary dM . This imposes the boundary condition SAf^ = on dM, which for purely magnetic 
solutions fixes the value of magnetic charge. 

The globally regular solutions have an arbitrary periodicity /3 of the Euclidean time coordinate. In the 
black hole case, the value of P is found by demanding regularity of the Euclideanized manifold as r ^ rh- 
It can easily be verified that the Hawking temperature expression Th = 1/(3 found in this way coincides 
with that given by the surface gravity computation. 

The physical interpretation of this formalism is that the class of re gula r stationary metrics forms an 
ensemble of thermodynamic systems at equilibrium temperature Th jSg. Z has the interpretation of 
partition function and we can define the free energy of the system F = — ^ogZ. 

Therefore 

log Z = -I3F = S~ PM, (49) 

or 

S^pM^hu (50) 

straightforwardly follows, with S the entropy of the system. 

To compute Icu we make use of the Einstein equations, replacing the R — 2A volume term with 2R\ 
WirCTf. For a purely magnetic ansatz {At = 0), the term T/ exactly cancels the matter field Lagrangian 
in the bulk action L,„ = —l/2Tr{Ffj^^F^'^). The Ricci component Rl is computed by integrating the 
Killing identity VVhXj, ~ RhcK'^, for the Killing vector if = 6f. The divergent contribution given by 
the surface integral term at infinity in R^ is also canceled by /surface + Icti yielding a finite expression of 
the action. For the metric ansatz (|24|l describing spherically symmetric and topological black holes static 
configurations one finds 

while the corresponding expression for static axially symmetric configurations described by the metric 
ansatz (fTH|) is 

The correponding expressions in the globally regular case are found by taking — s- in the above relations. 

Replacing now in H5()(l (where M is the mass-energy computed in Section 3.2, 3.3), we find S* = in 
the absence of an event horizon, while the entropy of the black hole solutions is one quarter of the event 
horizon expected. 
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E. The thermodynamics of n = 1 static black holes 



Based on the numerical results presented in Section III, we attempt here a discussion of the thermo- 
dynamic properties the A = — 3 EYM black solutions. To our knowledge these have not been previously 
considered. 

For simplicity, we will restrict our considerations to n = 1 static configurations {i.e. spherically symmetric 
or topological black holes). Thus we shall analyze black hole thermodynamics in a canonical ensemble, 
holding the temperature T and the magnetic potential at the boundary at infinity {i.e. the magnetic 
charge) fixed. The associated potential is the Hclmholz free energy F. 

The response function whose sign determines the thermodynamic stability is the heat capacity 



Stability follows from C > given the fact that black holes radiate at higher temperatures when they are 
smaller. 

The behaviour of the specific heat C can be easily understood from the state equation S = S{T) at 
fixed magnetic charge Qm- In Figure 11 we plot these curves for several values of Qm for fc = 1,0 and — 1 
black holes with unit winding number. 

For k — 1, the usual Schwarzschild-AdS behaviour is reproduced: the curves first decrease toward a 
minimum, corresponding to the branch of small unstable black holes, then increase along the branch of 
large stable black holes. This is in strong contrasts with the behaviour of the Abelian RNAdS solutions, 
which present only one branch and approach the extremal limit as T ^ (see Figure 11a). Note that 
the k = 1 solutions with small values around zero of the gauge potential at infinity appear to present a 
complicated thermodynamic structure. However these solutions are unstable to small perturbations and 
we will not consider them here. 

As seen in Figure llb,c, the heat capacity (|53|l is always positive for hairy black holes with zero or 
negative curvature horizon, (this appears to be valid for any wq). As a result, the fc = 0, 1 topological 
black hole solutions are always thermodynamically locally stable. 

It is instructive to plot also the free energy F = I / j3 as a function of temperature and various values of 
the magnetic charge (Figure 12). One can see again that for fc = 1 there are always two black hole radii 
associated with each temperature, for any value of Qm- Correspondingly, the smaller branch is unstable 
having negative specific heat. However, the action of the k — 1 solutions becomes positive for some critical 
value of the event horizon radius, for any value of the nonabelian magnetic charge. 

In the vacuum case, this indicates the existence of a phase transition. When the free energy is negative, 
the Schwarzschild -AdS black hole phase is dominant over the thermal AdS background phase. When 
the free energy changes its sign, the Hawking-Page phase transition between the AdS black hole and the 
thermal AdS background takes place 53] . 

As seen in Figure 8, for k — 1 monopole configurations the globally regular solution has minimal 
energy in its asymptotic class of solutions (with a given Qm) and so is the thermal background. Thus a 
phase transition should exist between the large black hole solutions and the corresponding globally regular 
configurations. 

In the absence of matter fields, the action of fc = 0, — 1 black holes is always negative. Therefore these 
configurations can also be globally stable and there is no phase transition. The inclusion of YM fields 
changes this behaviour and the low temperature solutions have a positive action. 



Restricting to static solutions, we find from — lim,.^oo ^h^^, the following background metric upon 
which the dual field theory resides 




(53) 



F. On the boundary CFT 



{d9^ + fk{0)dip^) - dt^ 



(54) 
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For k = 1 this is describes a 2 + 1 dimensional Einstein universe; for fc = it is 2 + 1 flat space, while 
k — —I describes the three dimensional open static universe. 

It would be desirable to compute some quantities in this background and to compare the results with 
the bulk predictions. The main problem is that, even in the vacuum case, the AdS4/CFT3 correspondence 
is much less understood that the AdS5/CFT4 case. For example, in Klebanov and Polyakov proposed 
a duality between a theory of massless higher spin gauge fields in AdS4 spacetime on the one hand and 
the 0{N) vector model at large N in three dimensions on the other. However, although the details of 
the boundary CFT will depend on the details of the bulk supergravity theory, the generic properties are 
expected to be independent of the precise features of the theory. 

Here we should remark that as found in Section 3, the nonabelian matter field in the bulk does not 
approach asymptotically a pure gauge configuration. For static solutions, the boundary form of the unit 
winding number nonabelian potential is 

= ^{Luonde + [^hMlr^ + ^^r^ M0)dcp}, (55) 

with wq a real constant. The corresponding boundary gauge field expression for k — 1 and a winding 
number n > 1 is 

A^°'> = i(l - wo){T^d0 -nsinOTj^dip}, (56) 

where t^(0, tp), Tg{9, ip) are suitable combinations of the Pauli matrices, whose form in given in Appendix 
A. One can see that the winding number enters the boundary gauge field expression. 

From the AdS/CFT correspondence, we expect the nonabelian hairy black holes to be described by some 
thermal states in a dual theory formulated in a background given by H54|) . The spherically and axially 
symmetric solitons will correspond to zero-temperature states in the same theory. This CFT will interact 
with a background SU(2) field given by (|^ . if^ . 

As conjectured in Ref. the dual field theory is the field theory of a stack of N coincident M2 

branes with an background external field coupled to the R-symmetry current. The bosonic sector of the 
Lagrangian for a single M2 brane is i2j| 

1 = j Sx^/^{Da^D°-^ + ii?$2), (57) 

where D = d — lA*^") is the gauge covariant derivative, and R is the Ricci scalar of the boundary CFT 
metric ^ab- 

Computing quantum effects for a generic SU(2) background field is a difficult task. Simpler results are 
found only for the case wq = 0. For spherical symmetry, the expressions (|55|l . (|56|l describe the field of a 
Dirac monopole with n units of magnetic charge 

^(") = cos 6ld^. (58) 

(For 71 > 1 this is proven by transforming the gauge connection IjASp to a special gauge such that has 
only a T3 component.) 

A computation of the effective action If,ff for a charged singlet scalar field propagating in a zero- 
temperature Euclideanized d = 3 Einstein-universe background and interacting with the U(l) field given 
by H58|l is presented in Appendix E, based on a zeta-function approach. From the basic relation we 
expect leff to present a qualitative agreement with the corresponding bulk computation. First we note 
that a straightforward evaluation of HE3p . HE14|I give positive values for Ieff{n), as expected from the bulk 
results (with the bulk action of globally regular solutions Ib = (3M). It is also interesting to compute 
the ratio Ieff{n)/Ieff{i) and to compare with the bulk results. One finds e.g. /e//(2)//e//(l) ~ 3.89, 
/e//(3)//e//(l) ~ 8.52 while /b(2)//b(1) ~ 2.46, /s(3)//b(1) ~ 4.72. The discrepancy beteween the bulk 
and boundary results increases with n. 

However, on the AdS side, there are two distinct bulk configurations with zero temperature and magnetic 
charge n. First, there is the nonabelian soliton solution with charge n discussed in the previous Section. 
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The second solution is an extremal RNAdS black hole with zero temperature and the same value of 
magnetic charge. It is not clear how the CFT knows to distiguish between these different bulk solutions. 
In any case, this is not possible within the model H57|l . 

There is also another problem with the action principle H57|l . Since the d — 4 Einstein-Maxwell system 
has electric-magnetic duality 41], one expects this duality to also be manifest on the CFT side. However, 
its realization in d = 3 is a rather subtle question (see e.g. the discussion in 'i^). For example, the U(l) 



dual of ^ is a A^°^ = T3 c dt/2 with = n^J (-1 + Vl + l2n'^)/6 for an extremal RNAdS solution. We 

can also perform a computation similar to that in Appendix E, for an electric A^"^ . One can see that, 
in the zero temperature limit, the parameter c will not enter the final results, which are similar to the 
vacuum case. Thus we conclude that the model (|57f) is too simple to mimic the expected features of the 
boundary CFT. 

However, we can use the AdS/CFT "dictionary" to predict qualitative features of a quantum field theory 
in the background (I54|l . For example, the expectation value of the dual CFT stress-tensor can be calculated 
using the relation ^] 



hh'^'Tbc. (59) 

Applying this prescription to the n = 1 static solutions, we find the standard form for the stress tensor of 
a (2+1) dimensional CFT 

<r^'>=-^[3u'^u'' + r% (60) 

where it° = Sf. 

A similar computation can be done for an axially symmetric configurations in the bulk. Considering the 
more general rotating case, we find the field theory stress tensor 

< >^< + < >(™*) (61) 

where < >(''*) is a contribution which survives in the static limit 
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(62) 



where 




and — 9, — Lp, — t. Here M, /2 are continuous variables which encode the bulk parameters. 
< >(''°*) is the part of the CFT induced by the rotation in the bulk and has the expression 

< Tab >(™*^= C ( 1 ) , with C = — — (ji -f j2 sin^ 9) sin^ 9. (64) 

The CFT stress tensor is covariantly conserved and manifestly traceless. Even for static configurations in 
the bulk, a winding number n > 1 of the bulk configurations implies /2 7^ and thus a 0— dependence of 
the dual theory stress tensor (although the boundary metric is spherically symmetric). This is a unique 
property of AAdS gravitating nonabelian configurations, since the boundary stress tensor of an Abelian 
solution with the same global charges has /2 = 0. This also suggests the dual theory should also be 
sensitive to the integer n, and is much more complex than the simple model IjSTI) . 

The form (|61|l of the dual CFT stress-tensor is puzzling from yet another point of view, since there is no 
global charge associated with the parameter /2. However, the expression (|61|) means that the dual CFT 
is able to discern between n > 1 nonabelian and embedded Abelian bulk solutions with the same set of 
boundary data. 
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V. CONCLUSIONS 

In this paper we have discussed the basic features of the nonabeUan solutions of a EYM-SU(2) theory 
with a negative cosmological constant A — —ig^ (where g is the gauge couphng constant of the YM theory). 
This theory corresponds to a consistent truncation of d = 4 A/" = 4 S0(4) gauged supergravity. Except 
for the NUT-charged solutions, wc have considered the A ~ — version of all monopole configurations 
previously discussed in the literature for a generic value of the cosmological constant. This includes 
both particle- like globally regular and black hole solutions. These configurations are asymptotically AdS, 
possesing a regular origin or a regular event horizon. Apart from spherically symmetric solutions (or 
topological black holes) we presented arguments for the existence of A = — Sg^ static, axially symmetric 
solutions. These configurations have no counterparts in the abelian theory. They generalize to higher 
winding number the known spherically symmetric solutions, presenting a nontrivial angle-dependence of 
matter fields and metric functions. 

The main feature of the EYM AAdS configurations is the existence of a nonvanishing nonabelian mag- 
netic flux on the sphere at infinity. The thermodynamics of the static black hole solutions has been also 
discussed to some extent. Apart from static solutions, we discussed AAdS configurations with a nonvan- 
ishing angular momentum which have not been presented before in the literature. 

All known EYM asymptotically flat conflgurations are likely to present A = —?>g^ generalizations. Thus 
we expect the existence AAdS counterparts of the EYM configurations discussed in ji^l, satisfying a 
complicated angle-dependent set of boundary conditions at infinity. It would be interesting to construct 
AAdS nonabelian solutions which possess only discrete symmetries Ei^ and to find the corresponding 
boundary stress tensor. 

One should remark that all solutions discussed here may be uplifted to d = 11 supergravity. However, 
the solutions we discussed here are generically not supersymmetric. Supersymmetric solutions are likely 
to exist, but we expect them to present naked singularities (this is the case of the abelian counterparts 
with the same amount of symmetry). Note that the planar BPS solution of the A = —3.9^ EYM model 
found in closed form in ^12] has a naked singularity. 

Apart from a discussion of the physical properties of various bulk EYM configurations, we attempted 
a preliminary discussion of these solutions in an AdS/CFT context. The A — —3g^ EYM configurations 
should give us information about the structure of a dual CFT in a background SU(2) field. However, a 
naive computation of the effective action of a charged scalar field in a U(l) magnetic background field 
gave some inconclusive results. Further progress in this direction would require a better knowledge of the 
structure of the dual theory defined on the 2-1-1 dimensional boundary metric. The results we discussed in 
the Section IV appear to indicate that the dual CFT is able to distinguish between various bulk solutions 
with the same set of boundary data at infinity. 

We think that these issues deserve further study, which should lead to a deeper understanding of the 
AdS/CFT correspondence. 
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APPENDIX A: THE AXIALLY SYMMETRIC GAUGE FIELD ANSATZ 

The construction of an axially symmetric YM ansatz has been discussed by many authors starting with 
Manton and Rebbi and Rossi [l^. The most general axially symmetric YM-SU(2) ansatz contains nine 
magnetic and three electric potentials and can be easily obtained in cylindrical coordinates a;^ = (p, (p, z) 

A, = \a^^\p,z)t; + \A^;f\p.z)T- + \a(^\p,z)t:, (Al) 
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where the only (/j-dependent terms are the SU(2) matrices (composed of the standard (ti, T2, T3) Pauli 
matrices) 



cos Ti + sm rup T2 , 



- suinip Ti + cos nip T2, 



T3- 



(A2) 



This ansatz contains an integer n, representing the winding number with respect to the azimuthal angle 
(p. 

Transforming to spherical coordinates r, ^, with (p = rsin^, z = rcos9), it proves convenient to 
introduce, without any loss of generality, a new SU(2) basis (t",t^,Tj^), with 



cos f 



cos ( 



The general expression (jAip takes the following form in spherical coordinates 



(A3) 



(A4) 



where A'j^dx'^ — A1:dr + Agdd + A^dip + Afdt. This ansatz is axially symmetric in the sense that a rotation 
around the z— axis can be compensated by a gauge rotation C^pA = D'^ 24], with \1/ being a Lie-algebra 
valued gauge function. For the ansatz (|A4|) . = ncosOr" /2 — nsm9Tg /2. Therefore we find F^^^ = -D^VF, 
where W ^ A^~'^. 

We use in this paper a reduced YM ansatz, employed also in all previous studies on EYM solutions, 
with five of the gauge potentials taken identically zero 



A 



(r) 



A 



A 



0. 



The consistency of this reduction can easily be proven at the level of the YM field equations. A suitable 
parametrization of the six nonzero components of which factorizes the trivial ^-dependence is |33) 



A^dx^" 



Hi 



dr + (1 - H2)de 



2 



nsmf 



i/3Y + (l-if4)f 



{dtp dt) 



T T 



dt 



(A5) 



One may consistently take O = in this ansatz; however the inclusion of this metric function simplifies 
the set of boundary condition for the rotating configurations. For fl — = Hq = 0, the static axially 
symmetric ansatz used in previous studied on AAdS static EYM solutions is recovered; further, by taking 
n = 1, Hi = H2, = 0, H2 = H4 = w{r) one finds a spherically symmetric static ansatz which reduces to 
l(T^ after a suitable gauge transformation. 

The ansatz (|A5|) satisfies also some additional discrete symmetries ^3 (i^ particular the parity 
refiection symmetry) and it is also invariant under Abelian gauge transformations U 



(A6) 



U = exp (^2^^r(r, 

To fix this residual gauge degree of freedom we choose the usual gauge condition 

rdrHi ~ deH2 = 0. 
The non-vanishing components of the field strength tensor are given by 
1 



Fl 



F. 



r 

sin 



r 

sin 6 



[Hi,0 + rH2,r] , 
rHi,r + Hi{H3 + cote)] , 
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F's^ = - sine [H-i^e + ifscot^ + HiHi - 1] , 
F' = Sin0[i/4,e+cot0(i^4-if2)-if2if3], 



et 



et 



FL = — 



Fl 



rH5 r + HiHe - - sin 9{Hi{l - H^) - H3 + rHs r) - ^ r sin ffiJg 
r 

- rHe r - HiH^ + - smeiHiHs + (1 - H^) + rH^ r) - ^ r sinOil ~ HA 
r [ r ' 

H^.e - H2H6 + - sin0(F2(l - Hi) - cotOHs - H^.e) ~ — sin^iJg 
r r 

He e + H2H5 - - sm0{H2H3 + cot6l(l - H^) -H^e)-— sin6'(l - HA) 
r ' r 



sin 9 



H^Hi + He{H3 + cote) - - sin6l(cot6i(l - H^) + H3) 



(A7) 



One can easily verify the matter ansatz is compatible with the metric form (|21|l , since the energy momentum 
tensor satisfies Ttr = TtB= T^r =Tip0 = 0. 



APPENDIX B: BOUNDARY CONDITIONS 



1. Static axially symmetric solutions 



These solutions are obtained for a truncation = of the ansatz ljA5|l and a metric given by (|18|l . 

To obtain AAdS axially symmetric configurations with a regular origin or event horizon and with the 
proper symmetries, we must impose the appropriate boundary conditions. The boundaries of the system 
are the origin/event horizon and spacelike infinity, the 2-axis and, because of parity reflection symmetry 
satisfied by the matter fields, the p-axis. The boundary conditions at infinity and along the z— and the 
p-axis (i.e. for 6 — 0^ '^f'^) are similar for both globally regular and black hole solutions. 

We start by setting the boundary conditions at infinity compatible with the AAdS asumption. For the 
metric functions one imposes 

/|r— 00 — ^|r— 00 — ^|r— C30 — f ; (-^-^) 

while the boundary conditions for the matter part are 

H2\r=oo — Hilr^oo — LOQ, ^^l|r=oo — -ff3|r=oo = 0, (B2) 

where there are no obvious conditions on the value of luq- For a solution with parity reflection symmetry 
(the only type we consider in this paper), the boundary conditions along the axes are 

Hl\e=0,TT/2 = H3\g^Q^^/2 = 0, 9ei?2|e=0,7r/2 = (^6^^4161=0,^/2 = 0, (B3) 

def\e=o,7r/2 = 9em|e=o,7r/2 = del\e=o,Tr/2 — 0. 

Therefore we need to consider the solutions only in the region < 9 < Tr/2. Regularity on the z— axis 
requires also 

H2\e=o ^ H^le^Q, m|e=o = ^|e=o- (B4) 
For globally regular solutions, the boundary conditions imposed at the origin are 

drf\r=Q = 9,-m|r=0 = drl\r=Q — 0, i/2|r=0 — -ff4|r=0 = 1- Hi\r=o ~ i?3|r=0 = 0. (B5) 

The boundary conditions satisfied by the black hole solutions at the event horizon are 



!\r=rH — '^\r=rH — Ar=rh — 0> -ffl|r=rh — Oi drH2\r=rh — 9r^^3|r=rft — drH^r=rh — Oi 



(B6) 
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2. Rotating solutions 

These solution are found for metric form H21|l and the YM ansatz (jA5|l . Similar to the static case, we 
consider solutions with parity reflection symmetry i.e. < 9 < 7r/2. 

A systematic analysis reveals that the boundary conditions for the magnetic potentials Hi, ..,i?4 and 
the metric functions f, I, m presented in the static case remain valid in the presence of rotation. For the 
supplementary functions H^, Hq, Q we imposed 

^^|r=oo = 0, H5\r=oo =Vcos9, iJe |r=oo = siu 6* (B7) 

at infinity, and 

de^\0=Q^Tr/2 = 0, deH^le^o ^ Hz\g^^/2 ^ 0, i^6|e=o = 9ei?6|0=7r/2 = 0- (B8) 

on the axis. For globally regular solutions, the following set of boundary conditions is imposed at the 
origin 

H^lr^QsinO + H(i\r=QC0s9 = 0, 9r-f^5|r=o cos ~ ^r^^elr^o sin ^^ = 0, f2|,~o = 0. (B9) 
The black hole solutions are found imposing at the event horizon 

ThH^lr^r^ + COSeQh = 0, rhH(i\r=r^ - siu 6'f7h = 0, fl\r=r„ = ^h- (BIO) 



APPENDIX C: GENERAL RELATIONS 



Solutions of the field equations are also classified by the nonabelian electric and magnetic charges Qe 
and Qm- The definition of conserved currents and charges in a non-Abelian Yang-Mills theory is a problem 
approached by different authors in the last decades (see e.g. I48ll-j5l| '). with various solutions. A gauge 
invariant definition for the nonabelian charges was proposed in |5Clj| (see also H^) 

Qe^ ^ ded^\Fe^\, Qm ^ ^ j dedip\Fe^l (CI) 

where the vertical bars denote the Lie-algebra norm and the integrals are evaluated as r — )■ oo. The 
expression for the magnetic charge implied by this definition is Qm = \k — ujQ\Vk/4:TT for spherically 
symmetric and topological black hole solutions and Qm = "-11 ~ ^ol ^'^^ axially symmetric configurations. 

The magnetic charge defined in this way is equal (up to a sign) to the expression found by using the 
usual (gauge dependent) definition 

fdSkV^ Tr{F''' T}, (C2) 



47r _ 

(with T = Tj, for the gauge ansatz (|27|l and T — Tr for the axially symmetric generalization (^Sj). 
In evaluating the electric charge expression one uses the asymptotic expansion of the electric potential 

^ cose{V + {cisin^ 9 + C2)/r), Hq sin9{V + [cssin^ 9 + Ci)/r), (C3) 

(with Ci real constants). 

The energy density of the solutions is given by the tt-component of the energy momentum tensor T^. 
Of interest here is the electric part of this component, Tr{Ffj_tF'^'^'^} and its integral 



= j Tr{F^tF'''}^d\, (C4) 



which measures the contribution of the nonabelian electric field to the mass/energy of the system. Similar 
to the purely abelian part, by using the YM equations JSJl this integral can be expressed as a total divergence 

-Me= ( Tr{F^tF>''}^d'x^ I Tr{AtF'''}dS^ - I Tr{AtF'''}dS^. (C5) 

J -/ 00 J eh 
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Thus, for globally regular configurations, a vanishing magnitude of the electric potentials at infinity implies 
a purely magnetic solution. In contrast, one finds rotating black hole solutions with V — which are 
supported by the event horizon contribution. 

Here we remark that the angular momentum of any solution admits also an expression in terms of surface 
integrals of the matter fields ,2^ . This can easily be proven by using the existence of a potential W ^ 
and the YM equation 

J= f T^^cPx = j 2Tr{Fr^F''' + Fe^F^*}yf^(fx = (f 2Tr{WF^''}dSf, - (f 2Tr{WF'''}dS^C6) 

J J J oc J eh 

The above relation takes a particularly simple form for globally regular configurations. By using the gauge 
field asymptotics, the total angular momentum in terms of matter field coefficients is 

J= — ( --Cl-C2+t^0(gC3 + 2c4)j . (C7) 

In the black hole case, (|C7|) should be supplemeted with the event horizon contribution. This relation can 
also be used to test the accuracy of the numerical results. 



APPENDIX D: NUMERICAL ALGORITHM 



The spherically symmetric and topological black hole solutions can easily be found by using a standard 
differential equations solver. Starting with a suitable set of boundary conditions at the origin/event horizon 
the equations are integrated towards r — > oo using an automatic step procedure and accuracy 10~^^. The 
integration stops when the AdS spacetime asymptotics are reached with a prescribed accuracy. 

The situation is much more complicated for axially symmetric configurations, which requires the solving 
of seven (in the static case) or ten nonliniar partial differential equations. All axially symmetric solutions 
presented in this paper have been found by using a similar approach to that employed by Kleihaus and Kunz 
in their studies of AF nonabelian configurations. We use in the numerical procedure a suitable combination 
of the EYM equations such that the diferential equations for the metric and gauge functions are diagonal 
in the second derivatives with respect to r. These equations are then discretized on a (r, 9) grid with 
Nr X Ng points. The angular coordinate 9 runs from to 7r/2 and the radial coordinate goes from to 
some large enough value r^ax (typically Vmax — 2 x 10'^-;- 5 x 10'^). For any type of solution, we tested that 
the relevant quantities are insensitive to the cut off value Vmax- The grid spacing in the r— direction is non- 
uniform, while the values of the grid points in the angular direction are given by 6^ = (k— / {2{Ng — 1)). 
Typical grids have sizes 150 x 30 points. We monitored also the remaining Einstein equations which are 
not directly solved, assuring that they are satisfied with a reasonable accuracy. 

In this scheme, a new radial variable is introduced which maps the semi infinite region [rc,oo) to the 
closed region [0, 1] (with rc = or r/i). For the globally regular solutions, our choice for this transformation 
was X ^ r/{r + 1). For the derivatives this leads to the substitutions 

rF, — > x{l - x)F,, , r^F^r.r ((1 - xfF^,,, - 2(1 - (Dl) 

for any function F in the differential equations. For the black hole solutions, we employed a new coordinte 
X defined as a; = 1 — r/j /r, which leads to the following substitutions in the differential equations 

rF^r (1 - x)F^^ r^F^r^r ^ (1 - x)^F^,,, - 2(1 - x)F,^ (D2) 

for any function F. 

The resulting system is solved iteratively until convergence is achieved. All numerical calculations for 
axially symmetric configurations are performed by using the program FIDISOL (written in Fortran), based 
on the iterative Newton-Raphson method. A detailed presentation of the FIDISOL code is presented in 
|52|. This code requests the system of nonlinear partial differential equations to be written in the form 
P(r, ^, u, Ur, ug, Ur9, Urr, ugg) = 0, (whcrc u dcuotcs the unknown functions) subject to a set of boundary 
conditions on a rectangular domain. The user must deliver to FIDISOL the equations, the boundary 
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conditions, the Jacobian matrices for the equations and the boundary conditions, and some initial guess 
functions. The numerical procedure works as follows: for an approximate solution u^^\ P(u(^^) does not 
vanish. Next step is to consider an improved solution w*^^^ — u'^' + Am, supposing that + Au) — 0. 

The expansion in the small parameter Aw gives in the first order = + Au) « P(u^^^) + ^{u'^^^)Au . 

This equation can be used to determine the correction Au^^^ = Au. Repeating the calculations iteratively 
(^{(3) — u(3) ^^^(2) qIq^^ fj^Q approximate solutions will converge, provided the initial guess solution is close 
enough to the exact solution. The iteration stops after i steps if the Newton residual P{u'^^^) is smaller 
than a prescribed tolerance. Therefore it is essential to have a good first guess, to start the iteration 
procedure. Our strategy therefore is to use a known solution as guess and then vary some parameter to 
produce the next solution. 

To obtain axially symmetric solutions, we start always with the n = I solution as initial guess and 
increase the value of the relevant parameters slowly. For static solutions, the parameter we vary is the 
winding number n. The physical values of n are integers. Rotating configurations are found e.g. by 
increasing the magnitude of the electric potential at infinity (for black holes, we vary also flh)- The 
iteration is done in small steps and eventually converges with a good enough accuracy. Repeating the 
procedure one obtains in this way solutions for requested values of the relevant parameters. For some of 
the configurations, we interpolate the resulting configurations and use them as a starting guess on a finer 
grid. 

FIDISOL automatically provides also an error estimate for each function, which is the maximum of the 
discretization error divided by the maximum of the function 52]. For the solutions discussed in this paper, 
the typical numerical error for the functions is estimated to be on the order of 10^"^. The output of the 
code was analysed and visualised mainly with MATHEMATICA. 



APPENDIX E: CHARGED SCALAR FIELD IN A EINSTEIN UNIVERSE WITH A U(l) 

DIRAC MONOPOLE BACKGROUND 

We consider the following action principle for a nonminimally coupled scalar field interacting with a 
background U(l)-field 



= - J {D^(t>D^'(|)* + M^(l)(j)* + ^RU*) ^/^cPx, (El) 
where D = d — iA^^^ , Af is the scalar field mass and ^ determines the coupling with the scalar curvature 

R = 2k/e. 

The zeta function approach implies the computation of the the eigenfunctions (j)^ and the eigenvalues 
\m of the differential second-order selfadjoint operator A = —D'^'^Df^ + A'P + ^R. Thus we consider the 
series with s e C (the prime on the sum means that any possible null eigenvalues are omitted) 

as\A)^Y^\-^\ (E2) 

N 

As is well-known, this series converges provided Re s > D/2. It is possible to continue the above sum into 
a meromorphic function of s that is regular at s = |54| . In a path integral approach, the effective action 
for a scalar field can be formally expressed as the functional determinant of the operator A as 

Ieff^~\\nAet{A|^Ji% (E3) 

where /i is an arbitrary renormalization mass scale coming from the path-integral measure. This determi- 
nant however is a formally divergent quantity and needs to be regularized. 

In a zeta function renormalization framework, the regularized determinant reads 

lndet(A^2) = -Q'{Q\A)-C.{Q\A)\n^?. (E4) 

We note that since C(0| A) = in odd dimensions (which is our case), the dependence on the renormalization 
scale drops out. 
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The eigenvalue equation Acf)]^ = ^N't'N can be solved by using the ansatz 

0^ = e*("^-"*)i^(6'), (E5) 

where F{9) is a solution of the equation 

satisfying certain boundary conditions at the limits of the 6 interval. Here we will restrict to the k = 1 case 
of a zero-temperature 2 + 1-dimensional Einstein universe background and a U(l) field A^"' = Q„i cos 9dip. 
Note that m should be an integer for fc = 1. 

The substitution F ~ g/y/s'mO, 9 — 2x transforms (jE6|) into the quantum mechanical problem of the 
Schrodinger equation with a Poschl- Teller potential 

- S + - Qrr^r - + ((m + Q„0' - (E7) 

dx'' 4 sm X 4 cos^ a; 

= (4Q^ + 1 - 4^2(^2 ^ ^^2 ^ _ ^ 0^ 

whose solutions are well known, see e.g. [s^. This leads to an eigenvalue expression 

Ajv = + 1 (1 + \m + Q,n\ + \m- Q„| + 2n) ' + +^R-Ql^- 1/4, (E8) 
where {N — m,n,Lu). The corresponding eigenfunctions are 

(t)N{x) = (cos^)l™+Q"l+i/2(sin^)l™-'3'"lp(l™"'5"l'l™+'3"l)(cos6')e*('"'^-"*), (E9) 

with P'^'^''^\x) the Jacobi polynomials. 

After integrating over lu, we arrive at a sum on the form 

^ n— m— — 00 

Restricting to the case of a massless, conformally coupled field, the resulting zeta function can be written 



00 00 



^^''^-^Tis)^ kl ^ (n+^(l + |m + 0™| + |m-Q„|)r^-^'=+i (Ell) 

^ ^ k—O n— Om— — 00 

Neglecting an additional contribution arising as a result of interchanging the order of summation, one 
performs first the sum over (m, n), finding the simple approximate expression 

C(^) - V^T^ EZo (2Cif(2s + 2fc - 2, Q„, + f ) (E12) 

+(20„ + 1) {Ch{2s + 2fc - 1, Q„ + i) - Ch(2.s + 2fc - 1, + §)) 

with C^r(s,a) the Hurwitz zeta functions, which are meromorphic functions with a unique simple pole at 
s = 1. In deriving this relation we used also [ssf 

J2 [m + n + a]-' = C„{s - 1, a) - (a - 1)Ch(s, a). (E13) 

m,n— 
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The zeta function (|E12ll is analytic throughout the complex s— plane except for s = 3/2 — n (with n = 
0, 1, . . .) where simple poles appear, while C(0) = 0. The derivative of this function evaluated at s = is 



C'(o) 



r(fc- 



-2V^ 



Cff (-2, Qrn + I) + {2Q,n + 1) {(h{-1, Qrn + ^) - Ch{-1, + §)) 

2Cff (0, + i) + (2Q™ + 1) (*(Q™ + f ) - ^'(Q™ + i)) 



(E14) 



-)2k+i 



+ Sfc=0 (fc+2 

where ^'(a;) — d\ogT{x)/dx. 



2Cff (2fc + 2, Q„ + I) + (2Q„. + 1) {CH{2k + 3, Q,„ + i) - Ch(2A: + 3, Q,„ + f )) 
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Figure 2. 

The profiles of typical spherically symmetric globally regular monopolc solutions are plotted for several values of 
the parameter h which enters the expansion at the origin of the gauge potential. 
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Figure 3. 

Mass M is plotted as a function of uoq for globally regular gravitating monopole static solutions. The winding 
number n is also marked. The configurations with n = 2, 3 are ax;ially symmetric. 
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Figure 4a. 
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Figure 4b. 

The mass M and the angular momentum J of non-Abehan globally regular rotating solutions axe shown as a 
function on the parameter V (Figure 4a, for wo = 0.262) and the parameter wo (Figure 4b, V = 2). Also shown are 
the electric charge Qe and the contribution Me of the nonabelian electric field to the total energy of the system. 
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Figure 5. 

Typical profiels of n = 1 black hole monopole solutions. The solutions with k = 0,-1 correspond to topological 
black holes with nonabelian hair. 
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Figure 6a. 
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Figure 6b. 
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Figure 6c. 

The mass-paxameter M, the asymptotic value of the gauge function loq, the value of the metric function p(r) at 
the event horizon and the Hawking temperature are represented as a function of the value of the gauge function uj 
at the event horizon for typical monopole black hole solutions. 
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Figure 7b. 
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Figure 7c. The value of the iionabchan gauge potential at the event horizon wy^ is plotted as a function 
of the event horizon radius Vh for spherically symmetric and topological black hole monopole solutions and several 
values of the magnetic potential at infinity. 




Figure 8a. 
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Figure 8c. The black hole mass M is plotted as a function of event horizon radius for k = 1,0,-1 
black hole monopole solutions with n = 1 and several values of the magnetic potential at infinity. 
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Figure 9b. 

The mass M (Figure 9a) and the Hawking temperature and the entropy (Figure 9b) are plotted as a function of wo 
for static black hole monopole solutions with ru = 1- The winding number n is also marked. The configurations 
with n = 2, 3 are axially symmetric. 
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Figure 10. 

The mass M, the angular momentum J, the electric charge Qe and the contribution Me of the electric field to the 
total energy of the system of non-Abelian black hole rotating solutions are shown as a function on the parameter 
V for fixed values of rn, uio and 0,h- Also shown are the entropy and the Hawking temperature. 
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Figure lib. 
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Figure 11c. The entropy is plotted as a function of temperature for k = 1,0,-1 black hole monopole 
solutions and several values of the magnetic potential at infinity. Here and in Figure 12a, the k = 1 curve with 
Wo = corresponds to RNAdS abelian black holes with unit magnetic charge. 




Figure 12a. 
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Figure 12b. 
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Figure 12c. The free energy F = M — TS is plotted as a function of temperature for k = 1,0,-1 black 
hole monopole solutions and several values of the magnetic potential at infinity. 



